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ABSTRACT 

Context. The Orion A molecular complex is a nearby (420 pc), very well studied stellar nursery that is believed to contain examples 
of triggered star formation. 

Aims. As part of the Herschel Guaranteed Time Key Programme CHESS, we present the discovery of a diffuse gas component in the 
foreground of the intermediate-mass protostar OMC-2 FIR 4, located in the Orion A region. 

Methods. Making use of the full HIFI spectrum of OMC-2 FIR 4 obtained in CHESS, we detected several ground-state lines from 
OH + , H 2 + , HF, and CH + , all of them seen in absorption against the dust continuum emission of the protostar's envelope. We derived 
column densities for each species, as well as an upper limit to the column density of the undetected HjO + . In order to model and 
characterise the foreground cloud, we used the Meudon PDR code to run a homogeneous grid of models that spans a reasonable range 
of densities, visual extinctions, cosmic ray ionisation rates and far-ultraviolet (FUV) radiation fields, and studied the implications of 
adopting the Orion Nebula extinction properties instead of the standard interstellar medium ones. 

Results. The detected absorption lines peak at a velocity of 9 km s , which is blue-shifted by 2 km s with respect to the systemic 
velocity of OMC-2 FIR 4 (Vlsr = H-4 km s -1 ). The results of our modelling indicate that the foreground cloud is composed of 
predominantly neutral diffuse gas («h = 100 cnT 3 ) and is heavily irradiated by an external source of FUV that most likely arises from 
the nearby Trapezium OB association. The cloud is 6 pc thick and bears many similarities with the so-called C + interface between 
Orion-KL and the Trapezium cluster, 2 pc south of OMC-2 FIR 4. 

Conclusions. We conclude that the foreground cloud we detected is an extension of the C + interface seen in the direction of Orion KL, 
and interpret it to be the remains of the parental cloud of OMC-1, which extends from OMC-1 up to OMC-2. 

Key words. Astrochemistry - ISM: clouds - ISM: abundances - ISM: molecules - ISM: photon-dominated region (PDR) 



1. Introduction 

The Orion region contains the nearest giant molecular clouds to 
the Sun and has been the object of a great variety of studies in 
the past few decades. Within this region lies the Orion A molec- 
ular cloud, a very active site of star formation, at a distance of 

. -420 pc (Hirota et al. I20071 Menten et al. 120071 . Located in the 
interior of a 300 pc diameter superbubble of gas inflated over 
the past 10 Myr by various OB associations (Reynolds & Ogden 

■ II 9791 1. Orion A is exposed to external forces that are believed to 
compress it and trigger star formation (e.g. Shimajiri et al. 1201 11 1. 

The northern portion of the Orion A complex is associated 
with the emblematic Orion Nebula (M42) and the Trapezium 
OB stellar cluster. Immediately behind the nebula lies the Orion 
Molecular Cloud 1 (OMC-1), which contains the BN infrared 
(IR) cluster and the KL molecular clump (see also the sketch 
in Fig. U). An integral-shaped filament (ISF) of dense molec- 
ular gas extends towards the north and south of OMC-1, with 
a length of about 15 pc and a relatively small width of ~0.5 pc 
(Bally et al. 1 19871 Johnstone & Bally 119991 . A velocity gradient 
is observed along the filament starting from Vlsr = 4 km s in 
the southernmost edge, up to Vlsr = 12 km s _1 in the northern 
part, with an abrupt change in M42. This velocity gradient is in- 



terpreted by Kutner et al. ( 119771) as rotation of the entire cloud 
about an east-west axis located roughly at the position of OMC- 
1 . The northern region of the ISF comprises two main molecular 
clouds active in star formation, known as OMC-2 and OMC-3. 
They are bound to the north by the reflection nebula NGC 1977, 
and have been extensively investigated at several wavelengths 
(see Peterson & Megeath 2008 for a review). 

Here we report the discovery of a diffuse gas component in 
the foreground of OMC-2. The data are part of the Herschel 
Guaranteed Time Key Programme Chemical HErschel Surveys 
of Star forming regions (CHESS0; Ceccarelli et al. 1201 Oi l, which 
has conducted unbiased spectral surveys of eight star-forming 
regions, each representative of a particular aspect of the star for- 
mation process: evolutionary stage, mass of the forming star, 
and/or interaction with the surroundings. Among the CHESS 
targets there is OMC-2 FIR 4, an intermediate-mass protostar 
and the brightest sub-millimetre source in OMC-2 (Johnstone & 

1 Herschel is a European Space Agency (ESA) space observatory 
with science instruments provided by European-led principal inves- 
tigator consortia and with important participation from the National 
Aeronautics and Space Administration (NASA). 

2 http://www-laog.obs.ujf-grenoble.fr/heberges/chess/ 
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Bally 1 19991 Kama et al. in prep.). The rich spectral dataset pro- 
vided by the HIFI spectrograph has allowed us to detect several 
ground-state lines from OH + , H20 + , HF, and CH + , seen in ab- 
sorption against the dust continuum emission of the protostars's 
envelope. 

We describe these observations in Sect. [2] The spectra and 
derived physical quantities are presented in Sect. [3] Section|4]de- 
scribes the modelling of the observations. In Sect. [5] we discuss 
the possible interpretations of our results, and a final summary is 
presented in Sect. [6] 

2. Observations 

The observations towards OMC-2 FIR 4 (a 2 ooo = 05 h 35 m 26.97 s , 
5 2 ooo = -05 d 09 m 54.5 s ) were carried out with the HIFI spec- 
trograph on-board the Herschel satellite (de Graauw et al. 
120101 Pilbratt et al. 2010) in single-pointing, dual-beam switch 
mode and with a redundancy of 8. The wide band spectrom- 
eter (WBS), whose spectral resolution is 1.1 MHz, was used. 
The double sideband (DSB) spectra were reduced and decon- 
volved using HIPEn version 7.0, and were subsequently ex- 
ported to GILDASQ for analysis. The final spectra presented 
here are equally weighted averages of the H and V polarisations, 
smoothed to a velocity channel spacing of 0.5 km s . 

Table Q] lists the observational parameters for each absorp- 
tion line detected towards OMC-2 FIR 4: beam efficiency, beam 
width, and the measured lcr RMS in main beam temperature, 

T m b- 

Following the detection towards OMC-2 FIR 4, as part 
of a Herschel OT2 accepted project, we recently observed 
OH + (lo,i/2-Oi ; 3/ 2 ) towards OMC-3 MM6, located ~2 pc far- 
ther north, to evaluate the extent of the absorbing foreground 
layer. The observations were made on 22 February 2012 in 
DSB, single-pointing mode, and were subsequently analysed 
with HIPE version 8.1. The corresponding observational param- 
eters are also listed in Table Q] Even though the RMS level is 
lower than for OMC-2 FIR 4 and the continuum level is slightly 
higher, the line is not detected towards this source. 

3. Results 

3.1. Spectra 

Figure Q] shows the observed absorption lines towards OMC- 
2 FIR 4, all of them corresponding to ground-state transitions. 
Four OH + lines, three of which are blends of two different hyper- 
fine components, are detected with a signal-to-noise ratio greater 
than 3. On the other hand, o-H20 + (l i, 1,5/2-00,0,3/2) is detected at 
4cr, and H30 + , whose ground-state transition falls at 984.7 GHz, 
is undetected. Both HF(l-O) and CH + (l-0) display clear absorp- 
tion and reach the saturation level. 

Two features can be noticed from Fig. Q] (see also Table Q]). 
Firstly, the absorption lines peak around a velocity of 9 km s _1 , 
which is shifted by more than 2 km s from that of OMC- 
2 FIR 4 (Vlsr = H-4 km s _1 ). Secondly, similar studies re- 
cently performed towards other sources typically reveal absorp- 
tion over a wide range of velocities, often several tens of km s 

3 HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, which consists of the ESA, the NASA Herschel 
Science Center, and HIFI, the Photodetector Array Camera and 
Spectrometer, and the Spectral and Photometric Imaging Receiver con- 
sortia. 

4 http://www.iram.fr/IRAMFR/GILDAS 



(e.g. Gerin et al. l20T0l Neufeld et al. l20T0l Monje et al. 120111 
and in some cases they are even contaminated by blending emis- 
sion lines (e.g. Phillips et al. 1201 Oi l. This is not the case here, 
where the lines have relatively small measured line widths and 
apparently trace one single velocity component. 

The centre of the o-H.20 + (l 1,1,5/2-00,0,3/2) line is blue-shifted 
with respect to that of the OH + lines by 1 km s _I (see Table[TJ. 
This is most likely due to the uncertainties in the rest fre- 
quency measurements reported in the literature for H20 + (see, 
e.g. Table 1 in Ossenkopf et al. 1201 01) . Here we have taken the 
rest frequencies from Murtz et al. 1 19981 

As mentioned above, OH + is not detected towards OMC- 
3 MM6 down to the noise level of the observations. 

3.2. Column density estimates 

The molecular column densities of each species displaying ab- 
sorption lines were estimated assuming that the particles are in 
the ground state, since their collisional coefficients are unknown. 
Under this assumption, the column density, N mo i, is given by 

N mol = 8n—^- TdV LSR , (1) 

^ul gu C* J 

where A u \ is the spontaneous emission Einstein coefficient, g\ 
and g u are the statistical weights of the lower and upper energy 
levels, respectively, and c is the speed of light. 

Since the ground state of both OH + and o-H20 + are split 
into two levels with a negligible energy difference, we assume 
that the population of each is proportional to their corresponding 
level degeneracies, namely g\ = 4 for F — 3/2 and g\ — 2 for 
F = 1/2. We therefore included the corresponding factor in Eq.Q] 
to obtain the total column density of each of these two species. 

A sketch of the energy levels involved in the OH + N = — > 1 
transition is shown in Fig. [2] 

The optical depth at a given velocity channel, r, is given by 



\ ^cont / 

The resulting velocity-integrated optical depths and column 
densities are listed in Table [1] together with some parameters of 
the transitions. The table includes upper limits for the undetected 
H30 + transition. The errors of the derived column densities are 
of the order of 30% and include the contribution of the spectral 
noise around the lines as well as a 25% uncertainty on the data 
calibration (Roelfsema et al. l2012l l. 

In the case of OH + , we used the four lines that are clearly 
detected above 3<r, all of which yield very similar values. The 
final OH + column density we adopted is the arithmetic aver- 
age of the four column-density derivations, which amounts to 
/V*oh+ = (2.4 + 0.4) x 10 13 cirT 2 . We computed 3<t upper limits to 
the column densities of their corresponding species for the un- 
detected OH + and H20 + hyperfine components. They all exceed 
the values listed in Table [1] and are therefore consistent with 
them. 

The highest column density of the three O-bearing species 
corresponds to OH + , followed by H20 + , and finally H30 + , with 
Noh+/Nh 2 o+ = 3.1, and JVh 2 o+/Nh 3 o+ > 9 - 5 - 

4. Modelling 

To model the chemistry of OH + , H 2 0\ HF, CH + , and the un- 
detected H30 + , and to determine the physical properties that 
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Table 1. Parameters of the observed absorption lines. 



Line 


V 




'/beam 


n 

"beam 


ltrRMS 


T 

1 cont 




FWHM a 


J rdV LSR 


N mo] 




(MHz) 


(10~ 2 s-') 


(%) 


(") 


(mK) 


(K) 


(km s" 1 ) 


(km S" 1 ) 


(km s" 1 ) 


(10 12 cm" 2 ) 


OMC-2 FIR 4 


CH + (l-0) 


835137.5 


0.64 


74.7 


26 


41 


0.48 


9.5 ±0.1 


6.1 ±0.2 


12.0 ± 1.6 


34.4 ± 4.5 


OH + (l ,l/2-0i,3/2) 


909158.8 


1.05 


74.4 


24 


33 


0.54 


9.3 ±0.1 


2.6 ± 0.4 


1.1 ±0.4 


23.1 ±7.4 


OH + (l 2 ,5/2-0l,3/2) b 


971803.8 


1.82 


74.1 


22 


95 


0.58 


9.2 ±0.1 


3.3 ±0.3 


5.0 ± 1.6 


23.8 ±7.3 


OH+(l M/2 -0 u/2 ) b 


1032997.9 


1.41 


74.1 


24 


52 


0.63 


9.4 ±0.1 


2.9 ± 0.4 


1.1 ±0.4 


24.9 ± 7.9 


OH + (l 1 , 3/2 -0 1 , 3/2 ) b 


1033118.6 


1.76 


74.1 


24 


52 


0.63 


9.3 ±0.2 


3.2 ± 0.5 


2.6 ±0.6 


22.7 ±4.8 


o-H 2 + (l u -0 ,o) d 


1115204.0 


3.02 


63.7 


19 


72 


0.68 


8.2 ±0.3 


4.4 ± 0.8 


1.8 ± 0.4 


7.6 ± 1.9 


H 3 + (0 i-loo) c 


984708.7 


2.30 


74.1 


22 


52 


0.59 






< 0.8 


<0.8 


HF(l-O) 


1232476.3 


2.42 


63.7 


17 


130 


0.75 


10.0 ±0.1 


2.8 ±0.3 


4.8 ± 1.4 


12.3 ±3.4 


OMC-3 MM6 


OH + (l .i/ 2 -0 1 _3 /2 ) b 


909158.8 


1.05 


74.4 


24 


13 


0.58 






< 0.1 


< 3.0 



a Measured from a Gaussian fit to the line 

b Blended with another hyperfine component (see Fig.flJ; A' rao i estimated assuming that all flux comes from this component 

c Undetected: 3o~ upper limits are given for JrdVLSR and A'moi, assuming FWHM = 2.6 km s _I 

d F = 5/2 - 3/2 transition; rest frequency from Miirtz et al. l 19981 
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Fig. 1. Spectral lines detected in absorption towards OMC-2 FIR 4. The green vertical lines mark Vlsr = 9.3 km s , i.e. the 
systemic velocity of the OH + lines. 
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Fig. 2. OH + energy level digram for the detected transitions. 
Solid arrows depict lines we observe at an S/N > 3, while dashed 
arrows mark the lines with an S/N < 3. 



best match our observations, we made use of the Meudon photo- 
dissociation (PDR) codeQ, developed by Le Petit et al. (2006). 
The model computes the chemical and thermal structure of a ho- 
mogenous plane-parallel slab of gas and dust illuminated by a 
radiation field coming from one or both sides of the cloud. We 
adopted the standard carbon and oxygen chemical network, plus 
the network for the fluorine chemistry following Neufeld et al. 
(2009). The adopted elemental abundances are 1.3 x 10~ 4 (car- 
bon), 3.2 X 10~ 4 (oxygen) and 1.8 x 10~ 8 (fluorine). We consid- 
ered a slab illuminated by the interstellar standard radiation field 
(ISRF; 1.7 G , where G = 5.9 x 10" 14 erg cm 3 is Habing's 
ISRF) on one side -since we do not expect multiple UV sources 
to lie close-by, which might irradiate heavily on both sides- and 
by a stronger field on the other side. 

As can be seen in Fig. [TJ the CH + (l-0) and HF(l-O) lines 
are saturated and have larger widths, suggesting the presence of 
more than one kinematical component. For this reason we based 
the comparison between models and observations (see below) 
on the lines produced by OFF and H20 + , which are optically 
thinner and narrower and appear to trace a single velocity com- 
ponent. The upper limit on the H30 + column density also places 
a strong constraint on the best-fitting model. 

Our grid of models spans a sufficiently large parameter space 
to find the solution that best fits the observations, as described 
below. 



4.1. Grid of models 

We ran two grids of models, a first one (coarse grid) sampling 
a large parameter space in gas density («h), visual extinction 
(Ay), incident far-UV (FUV) radiation field on one side of the 
slab (hfuv)> an d cosmic ray ionisation rate and then a finer 
grid around the best fit from the coarse grid. 

We probed FUV field fluxes up to 1700 Go because strong 
radiation fields are expected in this region and are also consis- 
tent with previous work, which claims that large amounts of UV 
radiation are necessary to explain CO and C + observations to- 



wards the OMC-2 region (Herrmann et al. 1997 Crimier et al. 
2009 ). Similarly, we covered a large range of ( CI values, as pre- 
vious work suggested enhancement with respect to the canonical 
value (Gupta et al. 2010 ). The range of gas densities chosen cov- 
ers typical diffuse interstellar cloud values, and we also selected 
relatively low values of Ay because the slab does not display a 
distinct CO line component (Shimajiri et al. 201 1). The range of 
parameters of the two grids are listed in Table |2l together with 
other parameters, such as elemental abundances and grain prop- 
erties, which we kept fixed for all models. 

Table 2. Parameters of the Meudon PDR models 



First grid Second grid 


«h (cm 1 ) 


10, 100, 1000 100, 200, 300 


A v (mag) 


0.5,1,3 1,1.5 


m fuv (Go) 


1.7, 17, 170, 1700 170, 340, 510, 850, 1190, 1700 


<rcr(10- 17 S- 1 ) 


1, 10, 100,1000 10, 30 


Elemental abundances wrt H nuclei 


C 


1.32xl0- 4 


O 


3. 19x10^ 


N 


7.50xl0~ 5 


F 


1.80xl0~ 8 


S 


1.86xl0" 5 


Other fixed parameters 


Rw 


3.1 




1.9 x 10 21 mag cnT 2 


Initial T gRS 


70 K 


Grain radii 


3 x 10~ 7 - 3 x 10" 5 cm 



Cardelli et al. (1989) reported R v = 5.5 and N H /A V = 
2.9 x 10 21 mag cirT 2 towards the Trapezium stars within the 
Orion Nebula. Since we are not certain that these values can be 
extended to the OMC-2 cloud, we adopted standard interstel- 
lar medium values (see Table |2}. However, to see the effect of 
changing these two parameters on the results, we ran a number 
of test models around the best fitting one with the Trapezium Ry 
and Nn/Ay values. The differences are discussed in Sect. 14.41 

4.2. Best-fit model 

Considering the coarse grid of models, only two solutions are 
found to lie close to the observed column densities of OH + , 
H 2 + and H 3 0+: 



• Ay = 1 mag, «h = 10 cm 3 , g a 

• Ay = 1 mag, «h = 100 cirT 
1700 G . 



10" 17 s"\ kfuv = 170 Go 
<T cr = 10- 16 s~\ M FUV = 



http://pdr.obspm.fr/PDRcode.html 



The latter is much more likely, since the former implies hav- 
ing a 60-parsec-deep slab of material between OMC-2 FIR 4 
and us but with a remarkably low velocity dispersion. We there- 
fore discarded this and consider only the n — 100 cirT 3 solution 
below. A posteriori we also show that the the high FUV field 
irradiation implied by this solution is also consistent (Sect. 15721 . 

To better constrain the parameters of the absorbing slab, we 
then ran a second grid of models that covered a finer-tuned set 
of values around those of the best solution, in total 72 additional 
models (Table|2] right column). The models with Ay = 1.5 mag 
and/or «h > 100 cm -3 predict a higher amount of H30 + than the 
upper limit we measured. On the other hand, the Ay = 0.5 mag 
models from the first coarse grid underpredict the observed 
amount of OH + . This leads us to conclude that the visual ex- 
tinction, Ay, of the foreground slab lies just around 1 magnitude 
and is not denser than 100 cm 3 . 
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The £cr = 10" 16 s" 1 models predict lower column densities 
than those observed, whereas the = 3 x 1CT 16 s" 1 models 
match them within the errors. 

In Fig. |3]we present a plot of the ratio of observed to mod- 
elled column densities, A^ bs/-^Vmodei> as a function of mfuv for a 
100 cm 3 dense slab for which £ cr = 3 X 10" 16 s" 1 . While OH + 
and H20 + converge around N bs /Wmodei = 1 f° r the whole range 
of mfuv plotted, the models overpredict the amount of H30 + for 
values of mfuv below ~ 1200 Go- We are therefore looking at a 
diffuse cloud of gas that is being heavily irradiated by an external 
FUV field. 



Table 3. Column densities predicted by the best-fit model with 
«fuv = 1700 Go (see text) 




100 



1000 



Ufuv (1.7 Go) 



Fig. 3. Ratio of observed to modelled column densities, 
AUs/Wmodei, of OH + (red), H 2 + (green) and H 3 + (blue), as a 
function of the FUV field for a slab with Ay = 1,«h = 100 cm" 3 , 
and £ cr = 3 X 10~ 16 s" 1 . The widths of the red and green bands 
represent the errors around the measured column densities of 
OFF and H20 + . The blue area marks the only range of FUV 
field values possible for the models to be compatible with the 
derived 3cr upper limit to the H30 + column density. A horizon- 
tal black dotted line is depicted at iVobs/A'modei = !■ 



Table [3] lists the column densities predicted by the best-fit 
model (mfuv = 1700 Go) for the observed species and for other 
species of potential interest. The CH + column density predicted 
by the best-fit model is lower than the observed value by a factor 
~250, a remarkably high discrepancy. Indeed, none of the mod- 
els we ran (Table |2]i reaches as close as a factor 20 below the 
observed column density of CH + . This molecule is discussed in 
more detail in Sect. 15. II 



4.3. Two HF velocity components 

What appears to be the best fit for the O-bearing species fails to 
reproduce the amount of HF observed, which is ~20 times higher 
than the model prediction. However, a closer look at the HF( 1-0) 
spectrum in Fig.Q~]suggests the presence of two different velocity 
components: a weak component around Vlsr of the foreground 
slab, and a stronger component at 10.5 km s _1 , much closer to 
the systemic velocity of OMC-2 FIR 4. 

According to our modelling, the amount of HF in the slab 
is around 10 12 cm" 2 . The remaining 10 13 cm" 2 (see Table [TJ 
presumably come from the denser OMC-2 FIR 4 boundary. 
Indeed, we ran a few test PDR models for higher density gas 
(riu = 10 4 cm" 3 ) and high FUV irradiation fields (between 170 



Species 






TVmoi (cm l ) 






observed 


modelled (ISM) 


modelled (Orion) 


H 2 




— 


7 — ■ , „ 1(1 

6.1 x 10' ' 


3.1 x 10 19 


H 3 + 




— 


1.5 x 10 12 


2.4 x 10 12 


OH + 


(2.4 


±0.4) x 10 13 


2.9 x 10 13 


6.1 x 10 13 


H 2 + 


(7.6 


± 1.9) x 10 12 


7.5 x 10 12 


7.2 x 10 12 


H 3 + 


< 


8.0 x 10" 


3.5 x 10" 


1.2 x 10" 


C+ 




— 


2.5 x 10 17 


3.8 x 10 17 


C 




— 


1.0 x 10 13 


7.0 x 10 12 


CO + 




— 


7.5 x 10 7 


6.1 x 10 7 


CO 




— 


3.1 x 10 9 


7.4 x 10 8 


HCO + 






1.4 x 10 s 


3.7 x 10 7 


CH 






4.7 x 10 s 


2.6 x 10 8 


CH + 


(3.4 


±4.5) x 10 13 


7.7 x 10 10 


1.3 x 10" 


C 2 H 






3.0 x 10 3 


2.8 x 10 2 


N 2 H + 






1.1 x 10 5 


1.3 x 10 5 


CN 






2.4 x 10 6 


1.2 x 10 6 


HCN 






2.7 x 10 3 


1.3 x 10 3 


HNC 






4.0 x 10 2 


6.6 x 10 1 


HF 


(1.2 


± 3.4) x 10 13 


7.1 x 10" 


2.9 x 10" 



and 1700 Go) and found column densities of 2 x 10 13 cm" 2 for 
HF. This is slightly higher than the value we measured, which 
may be explained by the fact that HF is depleted onto dust grains 
in denser molecular gas, as found by, e.g., Phillips et al. ( 120101 1. 
As for the O-bearing species, our test models predict low col- 
umn densities, of ~ 10" cm" 2 , consistent with having most of 
the flux of the observed lines coming from the slab and not from 
this red-shifted component. 

To summarise, our observations are compatible with the 
presence of two components: the foreground diffuse gas slab at 
a velocity of ~9 km s" 1 , responsible for most of the OH + and 
H20 + absorption, and the denser boundary of the OMC-2 FIR 4 
core, where most of the observed HF is produced. 



4.4. Reliability of the modelling 

The reliability of the best-fit model is subject to the uncertainty 
in the values of some input parameters. In particular, as men- 
tioned in Sect. 14. II we investigated how the extinction properties 
affect the predicted column densities by running a number of test 
models around the best-fit values. 

The last column in Table [3] lists the column densities of sev- 
eral species according to our best-fit model when adopting the 
values of Ry and Nn/Ay reported in Cardelli et al. ( 119891 1. These 
imply a deeper penetration of external dissociating and ionising 
photons and therefore a smaller amount of H2 in the cloud, which 
alters the formation of other molecules. For the species consid- 
ered in this study, this translates into an amount of HF lower 
by a factor 2.4, a higher amount of CH + by a factor 1.7, and a 
larger difference in the relative column densities of OH + , H20 + 
and H30 + , which causes OH + to be twice as abundant and H30 + 
three times less abundant. The column density of H20 + does not 
change significantly. 

In summary, the use of the Orion Nebula extinction proper- 
ties instead of the standard interstellar medium properties results 
in variations of at most a factor 3 in the predicted column den- 
sities. We find that, in order to better reproduce the observed 
column densities with the Ry and Nn/Ay values from Cardelli 
et al. (I1989K it suffices to adopt a slightly lower cosmic ray ion- 
isation rate of about 2 x 10" 16 s" 1 . The low density and visual 
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extinction, as well as the strong external FUV radiation field, 
need to be around the best-fit model values quoted in Sect. 14.21 

5. Discussion 

The results of our observations and modelling indicate that the 
OH + and H20 + lines, as well as part of the HF, are tracing a 
foreground cloud of diffuse gas (Ay = 1 mag, «h = 100 cm -3 ), 
exposed to a cosmic ray ionisation rate ( a — 3 x 1CT 16 s _1 and 
a strong external FUV field of 1200-1700 Go- This causes the 
heavily irradiated side of the slab to have a considerably higher 
gas temperature than the other side (Fig. [4]>. With this heavy ir- 
radiation, UV pumping might have an effect on the level popula- 
tions of the species under study, but given the lack of information 
about their collisional coefficients, we are not able to quantify it. 




\ 1 

0.0 0.2 0.4 0.6 0.8 1.0 
Av 

Fig. 4. Gas temperature as a function of visual extinction, Ay, 
predicted by the Meudon PDR code for a slab of gas with «h = 
100 crrT 3 , = 3 x 10~ 16 s _1 and irradiated from the right by a 
FUV field of 1700 G . 



Below we discuss the chemistry and geometry of the cloud 
and propose an interpretation to explain its nature. 

5.1. Chemistry of the foreground slab 

We illustrate here the effect of a strong external FUV field such 
as the one implied by our modelling and observations, on the 
chemistry and the resulting molecular abundances of the species 
studied in this work. Fig.|5]shows these as a function of Ay, as 
predicted by the Meudon PDR code for our best-fit model. The 
slab is exposed to an external FUV field of 1700 Go on one side 
(right), and to a typical FUV field of 1.7 Go on the other side 
(left). The predicted abundances are discussed in the following 
sub-sections, separated into families of species. 

5.1.1. Oxygen-bearing species: OH + , H 2 + , and H 3 + 

The strong incident FUV field affects the three O-bearing species 
in different ways: while it enhances the abundance of OH + at a 
certain depth, the amount of H20 + , and especially that of H30 + , 
are considerably reduced. 

The formation of these molecules is believed to begin with 
the production of H + and H3 + via cosmic ray or X-ray ionisa- 
tion of hydrogen, followed by charge transfer to produce + . 




0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 5. Abundance of OFF, H 2 + , H 3 + , CH + , and HF as a 
function of visual extinction, Ay, predicted by the Meudon PDR 
code for a slab of gas with «h = 100 cm 4 , ( a — 3 x 10~ 16 s _1 
and irradiated from the right by an FUV field of 1700 Go. 



Rapid hydrogen abstraction reactions of + with H2 follow to 
yield OH + , then H20 + , and end with the production of H30 + . In 
diffuse molecular clouds, which have high electron abundances, 
H30 + is destroyed via dissociative recombination to yield OH 
and H z O. 

This is our case: the heavy UV irradiation on the slab surface 
gives rise to a high abundance of protons that favour the produc- 
tion of OH + , but at the same time it destroys molecular hydro- 
gen, thus hindering the cascade of reactions with H 2 that leads 
to the formation of H 2 + and subsequently to H?0 + . The strong 
FUV field is therefore responsible for the relative column den- 
sities (A^oh+ > Nh^o+ > Nh 3 o+) we measure in the foreground 
slab. 

5.1.2. Hydrogen fluoride: HF 

The abundance of HF mirrors that of molecular hydrogen, with 
an estimated HF/H 2 abundance ratio within an order of magni- 
tude from 3.6 x 10~ 8 (Neufeld et al. |2005l Indeed, HF has a 
higher binding energy than that of H 2 and, as a result, fluorine 
atoms are unique in reacting exothermically with the dominant 
constituent of interstellar molecular clouds, H 2 , producing hy- 
drogen fluoride via the reaction 

F + H 2 -> HF + H 



HF is therefore predicted to be the main reservoir of gas- 
phase fluorine in diffuse clouds where the amount of molecular 
hydrogen exceeds that of atomic hydrogen. This is reflected in 
Fig- El where HF is less abundant close to the heavily irradiated 
surface of the slab, where the incident FUV photons dissociate 
H 2 , thus preventing the above reaction to take place and form 
HF. 

5.1.3. Methylidyne cation: CH + 

Unlike the O-bearing species and HF, the abundance of CH + 
does not change significantly with increasing depth within the 
slab (Fig. |5]l even though previous works (e.g. Stauber et al. 
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2004) predicted the abundance of CH + to be strongly enhanced 
with increasing UV field. 

Our best-fit model for the foreground cloud (Sect. 14.21 ) pre- 
dicts an abundance of about 2 x 10" 11 with respect to the total 
amount of hydrogen. As mentioned in Sect. 14.21 this is lower by 
a factor ~250 than the observed value, a remarkably high dis- 
crepancy. It is likely that, as in the case of HF, there are several 
kinematical components hidden in the CH + (l-0) line, although 
given its level of saturation, these are impossible to separate. 

Nevertheless, if we assume that the entire CH + absorption 
comes from the foreground slab and adopt A^h = 1.9x 10 21 cm" 2 , 
which corresponds to that of the best-fit model, the CH + mea- 
sured abundance is as large as 1.8 x 10" 8 . This value is of the or- 
der of that found by other authors in analogous studies (Bruderer 
et al. [20101 Falgarone et al. I2010al 1201 Obi ) , demo nstrating, as 
stated also for instance by Falgarone et al. (2010a) and Sheffer 
et al. (120081) . that the CH + abundances observed in the diffuse 
interstellar medium are several orders of magnitude higher than 
the predictions of UV-driven steady-state models, like the one 
we have used here. Non-equilibrium chemistry models need to 
be employed to correctly predict the CH + column densities. For 
example, Federman et al. ( 119961) and Sheffer et al. (f2008 ) repro- 
duced the observed CH + abundances by including non-thermal, 
magneto-hydrodynamic motions in their models, while Godard 
et al. (120091 ) achieved similar results by including supersonic tur- 
bulence in their chemistry modelling of diffuse interstellar gas, 
and illustrated in their Fig. 8 how standard PDR models fail to 
match the observed values, underestimating them by orders of 
magnitude. 

Finally, it is important to bear in mind that the failure of the 
Meudon code to predict CH + column densities has effects on 
other families of species. For O-bearing molecules, the forma- 
tion of CH + leads to CO via the following reactions (Federman 
et al. [19961 ): 

CH + + O CO + + H 

CO + + H ^ CO + H + 

Subsequent charge exchange between H + and O will lead to 
an enhanced amount of + , thus affecting the oxygen chemistry 
described in Sect. 15.1.11 The inability of standard PDR models 
to correctly reproduce CH + abundances is therefore an impor- 
tant caveat for using such models, since it has deeper implica- 
tions on the chemistry of other species, in particular the oxygen 
chemistry through the above-mentioned channel. 

5.2. Geometry of the slab 

An approximate calculation allows us to check whether the high 
FUV field necessary to explain our observations (between 1200 
and 1700 Go) is consistent. We assume the nearby Trapezium 
OB association, at a projected distance of 2 pc, to be responsible 
for this heavy irradiation. The Trapezium FUV output is domi- 
nated by four young stars: 8 Ori C, with a spectral type between 
04 and 06 (Mason et al. [19981 Baldwin et al. [19911 ). and three 
early B-type stars, 9 Ori A, B, and D. Table |4]presents the effec- 
tive temperature (from Panagia fl973l ) and luminosity integrated 
over the FUV range, namely between 912 and 2400 A, for each 
of the four stars. 

A few parsec north of OMC-2 FIR 4 lies the NGC 1977 Hn 
region, where the earliest-type ionising star is a Bl, which is 
very faint in FUV output relative to the Trapezium. We there- 
fore consider only the contribution of the four main stars of the 



Table 4. Properties of the dominant Trapezium stars and distance 
from the absorbing slab 



Star 


Spectral type 


Teff 


£fuv 






(10 3 K) 


(10 3 Lo) 


Ori A 


B0.5 


26.2 


13 


eOriB 


Bl 


22.6 


6.3 


eOriC 


06-04 


42-50 


230 - 1900 


(?OriD 


B0.5 


26.2 


13 




Distance (pc) 






"fuv = 1200 Go 


"FUV 


= 1700 G 


6 Ori C: 06 


2.0 




1.7 


G Ori C: 04 


5.4 




4.6 



Trapezium. Depending on whether 6 Ori C is an 06 or an 04 star, 
this results in a distance ranging from 2.0 to 5.4 pc between the 
Trapezium and the slab, for an irradiation of 1200 Go, and from 
1.7 to 4.6 pc for 1700 Go (Table|4](. This is a rough estimate be- 
cause we neglected the unknown contribution of dust extinction, 
which would yield a smaller distance. Despite this caveat, our 
calculation suffices to conclude that the high external FUV irra- 
diation implied by our modelling and observations is consistent 
with the projected distance of 2 pc between OMC-2 FIR 4 and 
the Trapezium. 

According to our best-fit model, the slab has a visual extinc- 
tion of 1 magnitude or, in other words, a total hydrogen column 
density of 1.9 x 10 21 cm" 2 . Taking into account its low density 
(«h = 100 cm" 3 ), this implies a line-of-sight size of 6 pc. 

5.3. Nature of the slab 

Numerous studies have been performed to investigate the fore- 
ground material towards OMC-1, which contains the Orion KL 
star-forming region. These are mostly based on UV data of 
atomic and/or ionised metals, and HI absorption line observa- 
tions at 21 cm. More recently, Herschel observations have pro- 
vided additional information (see e.g. Bergin et al. 120101 ). In con- 
trast, no analogous research has been carried out towards OMC- 
2 FIR 4 and the rest of the OMC2/3 region, mostly because of the 
lack of strong UV background sources like the Trapezium in this 
direction. We can nevertheless compare our results with those 
found in past studies for OMC-1 in an attempt to search for links 
and understand the nature of the foreground cloud we observe. 
The foreground material towards the Trapezium and OMC-1 can 
be grouped into the following two components: 

1 . Orio n's Veil: This term was used by O'Dell (120011 ) and Abel 
et al. (2004 2006) to denote a cloud of predominantly neu- 
tral gas that lies between the Trapezium cluster and us. It is 
also known as the "foreground lid" (e.g. Wilson et al. 119971 ). 
Orion's Veil has been observed and mapped in HI at 21 cm 
by several authors (e.g. van der Werf & Goss 1990, Wilson et 
al. 119971) . and neutral and ionised metals have also been de- 
tected in absorption at UV wavelengths (reviewed in O'Dell 
1200 II and Abel et al. 2006). This system appears to have two 
main components, named A and B, at 1 km s" 1 and 5 km s , 
respectively. Photoionisation models performed by Abel et 
al. (2004 2006) located the Veil at a distance of about 2 pc 
in front of the Trapezium, with a density of ~10 3 cm" 3 and a 
thickness of 1 pc. 

2. C + interface between Orion KL and the Trapezium: This 
component has been observed in carbon recombination lines 
tracing a C + region at 9 km s" 1 (Balick et al. 119741 ), and 
more recently in absorption at the same velocity in OH + and 
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Fig. 6. Sketch of the OMC- 1/2/3 region as it might have appeared before the formation of the Trapezium cluster (top), and how 
it looks at present time (bottom). The foreground cloud we detected is labelled "fossil cloud", and the colour gradient across it 
represents a density gradient roughly in the north-south direction (see text). The sizes of each component and their velocities 
relative to the systemic velocity of OMC- 1 and the foreground cloud (i.e. Vlsr = 9 km s -1 ) are marked. The observer, represented 
by a small picture of the Herschel satellite, is situated to the right of each panel. 



H z O + (Gupta et al. l20T0b . and HF (Phillips et al. EOTOb with 
Herschel-HIFI. The velocity is closer to the systemic veloc- 
ity of the molecular gas in Orion KL (Vi sl =8-9 km s _1 , 
eg. Shimajiri et al. 1201 11 1. This region is identified as a 
"C + interface" between the Orion KL molecular gas and the 
Trapezium Hn region. This component has strikingly similar 
characteristics to those seen towards OMC-2 FIR 4 in terms 
of detected species, velocity, and relative measured column 
densities. 

A summary of the properties of these foreground compo- 
nents, including OMC-2 FIR 4 from this work, is reported in 
Table E] 

The complexity of the velocity field in Orion A has been 
noted by several authors (see review by O'Dell 2001). In a sim- 
plified picture, it reduces to the following three Vlsr compo- 
nents: 



- Trapezium Hn region: -3 km s 

- Foreground lid: 1 and 5 km s _1 

- dense molecular gas (Orion KL): 9 km s _1 , with a south- 
north velocity gradient along the ISF (see Sect.[TJ. 

Wilson et al. (Q997) proposed a picture (see their Fig. 7) 
in which the Trapezium drives a strong wind, creating a cav- 
ity of low-density gas around it that moves towards us. This 
expansion is halted behind the Trapezium by the high-density 
molecular gas, but is allowed to extend in front of it, where there 
is no dense gas material to stop it. Beyond this cavity lies the 
ionised gas region, which is constrained between the foreground 
lid and the C + interface behind it. The foreground lid's receding 
velocity with respect to that of the Hn region is interpreted as 
a flow away from this neutral material towards the lower pres- 
sure ionised gas in which the Trapezium is contained. The C + 
interface, also observed in absorption against the background 
dust continuum emission of Orion KL in OH + , ILO + and HF, 
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Table 5. Properties of the foreground material towards OMC-2 FIR 4 and towards Orion KL. 





OMC-2 FIR 4 


Orion KL" 


Orion s Veil (A) 


Orion s Veil (B) D 


Species 


OH+, H 2 + 
HF, CH + , C + 


OH+, H 2 + 
HF, C + 


HI, metals 


HI, metals 


Vlsr 0™ s ~') 


9 


9 


1.3 


5.3 


«h (cm" 3 ) 


-100 


-1000 


300 


2500 


N ffl (cm- 2 ) 


-10 21 


-10 21 


1.6 x 10 21 


3.2 x 10 21 


Thickness (pc) 


6 


0.5 


1.3 


0.5 



<n) Gupta et al. (2010), Phillips et al. ( [20101 
(b) Abel et al. 120041120061 and references therein 



as mentioned above, is confined between the Trapezium and the 
molecular gas, which are about 0.5 pc away from each other 
(Zuckerman l 19731 1. The results found by Gupta et al. (2010) sug- 
gest a density of about 10 3 cirT 3 for this interface, and possibly 
a high ionisation rate (£ ~ 10~ 14 s~'). 

The similarity of our slab properties (velocity, column den- 
sity) to those of the C + interface component lends support to a 
scenario in which our slab is an extension of that neutral gas re- 
gion. The slab is blown-up to a lower density (100 crrr 3 ) and 
larger line-of-sight size (6 pc) by to the absence of elements that 
would confine it into a smaller and denser area, as for the C + 
interface, which is strongly affected by the Trapezium. 

The absence of OH + absorption towards OMC-3, after com- 
parison between the models and the derived upper limit on Now 
(Table [TJ, can be interpreted by either of the following two pos- 
sibilities: (i) the foreground cloud does not extend much far- 
ther to the north of OMC-2; (ii) it does reach the foreground 
of OMC-3 but is shallower (Ay < 0.5) and/or more weakly irra- 
diated (mfuv ~ 1, 200 Go) than the gas seen towards OMC-2. 
This is consistent with the larger distance of OMC-3 from the 
Trapezium cluster. 

The cosmic ray ionisation rate implied by our observations 
(£ = 3 x 10~ 16 s -1 ) is consistent with the typical value measured 
in diffuse clouds (Indriolo et al. 120071) . 

We stress that the foreground gas we detect towards OMC- 
2 FIR 4 has not been reported previously and therefore adds to 
the large ensemble of components present in the Orion A com- 
plex. We interpret this foreground gas seen towards both OMC- 
1 and OMC-2 as the remnants of the OMC-1 parental molec- 
ular cloud, which once contained the primordial material that 
later allowed the formation of the Trapezium, Orion KL, and the 
BN infrared cluster. Figure [6] presents two schematic cartoons 
of how the OMC- 1/2/3 complex might have looked like before 
the Trapezium formed {top) and how it currently appears accord- 
ing to our interpretation of the foreground slab, labelled "fossil 
cloud" (bottom). 



6. Summary 

As part of the Herschel key programme CHESS, we have anal- 
ysed a collection of absorption lines seen towards OMC-2 FIR 4. 
These include the ground-state transitions of OH + , H20 + , HF, 
and CH + . Our results can be summarised as follows: 

1. The detected absorption lines peak at 9 km s _1 , blue-shifted 
by 2 km s with respect to the systemic velocity of OMC- 
2 FIR 4 (Vlsr = 11-4 km s~'), indicating the presence of a 
foreground layer of gas in this direction. 

2. Analysis using the Meudon PDR code shows that, regard- 
less of whether we adopt standard interstellar medium ex- 
tinction properties or those reported in the literature for the 



Orion Nebula, the slab is composed of predominantly neu- 
tral diffuse gas («h = 100 crrr 3 ), which is heavily irradiated 
(~ 1500 Go) by an external source of FUV most likely aris- 
ing from the nearby Trapezium OB association. The slab is 
6 pc thick and bears many similarities with the so-called C + 
interface between Orion-KL, 2 pc south of OMC-2 FIR 4, 
and the Trapezium cluster. 
3. Combining this work and previous studies carried out to- 
wards Orion KL, we conclude that the foreground slab we 
discovered is the extension of the C + interface seen in the di- 
rection of Orion KL, and interpret it to be the remains of the 
parental cloud of OMC-1, which extends from OMC-1 up to 
OMC-2. 
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